A simple method for obtaining electron scattering phase shifts from energies of an 

atom in a cavity 



in 
o 
o 

(N 

>v 



(N 



I.M. Savukov 

Department of Physics, Princeton University, Princeton, New Jersey 08544 
(Dated: February 2, 2008) 

We present a simple method for obtaining elastic scattering phase shifts and cross sections from 
energies of atoms or ions in cavities. This method does not require calculations of wavefunctions 
of continuum states, is very general, and is extremely convenient from practical point of view: 
some conventional computer codes designed for the energies of bound states can be used without 
modifications. The application of the method is illustrated on an example of electron scattering from 
Kr and Ar. From Brueckner orbital energies in variable cavities, we have obtained ab initio cross 
sections that are in close agreement with experiment. The relativistic effects are also considered 
and found to be small below 10 eV. 

PACS numbers: 34.80.Bm, 31.15.Ar, 31.25.-v, 31.30.Jv 
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Conventional methods of calculations of scattering 
cross sections are cumbersome, inconvenient, and very 
often inaccurate. This is only because they all are based 
on computing continuum, or sometimes quasicontinuum, 
wavefunctions and asymptotic fittings to extract phase 
shifts. Such an approach requires modifications of con- 
ventional atomic structure codes, developed for bound 
states, or just writing new programs altogether. For a 
known potential it is not a difficult task - this is why 
numerous semi-empirical calculations can be found in 
the literature - but the level of accuracy and theoreti- 
cal uncertainty of calculations based on ad hoc poten- 
tials can not be totally satisfactory. For ab initio cal- 
culations already complicated codes have to be rewrit- 
ten, which takes considerable amount of time. For multi- 
configuration Hartree-Fock (MCHF) method this was un- 
dertaken by Saha 0,0 to obtain ab initio results in agree- 
ment with experiment. However, many-body perturba- 
tion theory (MBPT) methods, which were developed for 
fundamental symmetry tests, have not been used for cal- 
culations of electron scattering cross sections. 

The method we propose in this letter is very simple and 
general: instead of finding continuum wavefunctions and 
fitting them to asymptotical solutions to obtain phase 
shifts for given electron energies, we impose a boundary 
condition on an atom, an ion, or a molecule to make 
the spectrum discreet and then from discreet energies 
extract phase shifts which are uniquely related to these 
energies. Thus the problem of phase shifts is converted 
into a conventional problem of finding energies of bound 
states. Especially simple relation exists, as we will show, 
in the case of an atom in a spherical cavity. 

It can be shown that continuum and quasicontinuum 
wavefunctions are equivalent. For example, in Ref. it 
was stated that B-spline solutions obtained in a cavity 
can be interpreted as a representation of true continuum 
states with a different normalization, and the energy of 
the quasicontinuum states can be set to an arbitrary pos- 
itive value by adjusting the size of the cavity. There are 
also other methods that give B-spline continuum wave- 
functions at any energy: the Galerkin method least- 



squares approach 0, [(J , and free boundary condition ap- 
proach |2J. The emphasis in these works is placed on ap- 
plications of B-splines which are very often bundled with 
the cavity boundary conditions: for the method proposed 
here, however, the boundary conditions are more essen- 
tial than B-splines, which are still convenient for evalua- 
tion of radial integrals in high-precision MBPT calcula- 
tions p|. 

While our method can be justified mathematically in 
quite general assumptions, it is not yet obvious that the 
method will be accurate in practical calculations, so we 
will illustrate the usefulness and accuracy of the method 
on specific examples such as the MBPT calculations of 
electron elastic scattering phase shifts and cross sections 
from Ar and Kr. The MBPT is chosen because it can pro- 
vide the best accuracy for a negative mono-valent ion, e.g. 
Ar+e~, uses cavity-bound basis functions, so the code 
does not need modification, and allows systematic consid- 
eration of correlations. To obtain correct electron-noble 
gas scattering cross sections it is necessary to include the 
direct and exchange potentials from a frozen noble-gas 
atom as well as core polarization effects. The direct part 
of the unperturbed atomic potential produces a phase 
shift opposite to that of the exchange part. The combined 
phase shift from the frozen atom is opposite in sign to 
that of the core-polarization. We will show that so-called 
Brueckner-orbital (BO) approximation treats accurately 
these effects and results in good precision. The accuracy 
can be further improved by using all-order couple-cluster 
method, or other accurate methods developed for mono- 
valent atoms. The calculations of phase shifts from en- 
ergies for other systems should be also possible and will 
be undertaken in future. 

Apart from illustration purpose, the calculations will 
serve to provide accurate ab initio cross sections for com- 
parison with other theories and experiments and to im- 
prove understanding of this particular system. Despite 
longer than a century history [t| that experiments on 
electron interaction with gases have, many questions re- 
main open and this area of research is still very active. 
Elastic scattering of electrons on noble-gas atoms is of 
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particular interest since many precise measurements are 
available providing tests for theories which all with a few 
exceptions are not of ab initio type and are based on pseu- 
dopotentials to take into account exchange interaction 
and significant polarizability of noble-gas atoms by an 
electron. Although elaborate complicated semi-empirical 
effective potentials have been developed to achieve good 
accuracy of calculations, many different calculations and 
measurements are still in disagreement, and there is 
clearly significant uncertainty in theoretical understand- 
ing. This situation exist in almost all noble-gas atoms. 

For example, motivated by uncertainty in cross sec- 
tions at low energies, which are important for extraction 
of scattering lengths, first ab initio calculations of low- 
energy electron scattering from neon pj and argon Q 
based on MCHF method to account for polarization 
effects have been reported, and good agreement with 
experiment has been demonstrated. Although MCHF 
method is very effective in general for the considera- 
tion of complicated open-shell ions, in monovalent atoms 
and low-charge ions this method has lower accuracy than 
MBPT methods, which were not applied to calculations 
of electron scattering from noble-gas atoms, probably due 
to complications associated with continuum states. 

With the aid of partial wave expansion, 
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a total elastic cross section <r t can be found from phase 
shifts Si 



<Jt = pf X)CM + 1) sin 2 Si 



(2) 



which arc normally extracted from asymptotic behavior 
of radial wavefunctions Pi(r) obtained by numerical so- 
lution of radial Schrodinger equation 



d 2 Pi(r) 
dr 2 



U(r) 
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Pi{r) = (3) 



for a given energy E = k 2 /2 as a parameter. In this 
equation U(r) is some effective potential which describes 
approximately direct and exchange interaction as well as 
the attraction due to core polarizability. (Atomic units 
are used in all equations.) The radial wavefunctions can 
be also obtained by using ab intio atomic structure meth- 
ods such as MCHF 0, Q or MBPT. Because wavefunc- 
tions are not always available in precision MBPT calcu- 
lations and most codes output either energies or matrix 
elements, the extraction of phase shifts from wavefunc- 
tions is not very convenient. However, it is not necessary: 
phase shifts can be obtained from energies of an atom 
bound to a cavity, which is a natural setting in MBPT 
calculations. The extraction is possible because the cav- 
ity uniquely encodes phase shift information into energies 
of quasicontinuum states and quasicontinuum wavefunc- 
tions are proportional to true continuum wavefunctions 



if their energies are the same. The last statement can be 
easily proved since the continuum and quasicontinuum 
wavefunctions are both unique solutions of the radial dif- 
ferential equation with the same boundary condition at 
r — > 0, the same energies, although with different nor- 
malization conditions and maybe sign convention. The 
equivalence of quasicontinuum and continuum states was 
also stated in Ref. |3j. At large r continuum and quasi- 
continuum solutions approach asymptotically the solu- 
tion in empty cavity proportional to rji(r), where ji(r) 
are spherical Bessel functions, and the effect of the atomic 
potential is only in phase shifts which can be determined 
from the asymptotic form of the wavefunctions or from 
energies for a known cavity radius R, which is our pro- 
posed method: 
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where xi n is the n th zero of the spherical Bessel function 
jl (x) . Accurate values of xi n can be found in mathemat- 
ical reference books, for example on page 467 of Ref. ^3] , 
column j va where v = 1 + 1/2 and s = n. For 1 = 0, 
xon = tvk. The lowest quasicontinuum state of a given 
symmetry has to be used with the first zero of the cor- 
responding spherical Bessel function, the next state with 
the second zero, etc. 

Energies of quasicontinuum states are calculated in 
Brueckner-orbital approximation, which accounts for 
core-polarization effects with relatively high precision. 
First, the Dirack-Hartree-Fock (DHF) equation is solved 
for a closed-shell atom (Ar or Kr) . Then in the obtained 
DHF potential, B-spline finite basis is generated. In this 



basis, the Hamiltonian matrix hi 
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is calculated and diagonalized to obtain BO energies. 
The summation runs over core states c, excited states 
n, m, and angular momenta fc; the matrix elements are 
calculated between all possible states i and j. The cou- 
pled radial integrals X k (abcd) and Zk(abcd) are defined 
for example in The self-energy matrix elements 

Sjj(eo), which take into account dominant part of core- 
polarization effects, depend on electron energy £o and 
contains non-local interaction, so that they can not be 
approximated accurately with a single effective poten- 
tial unless the energy range is small, e < Sq, and the 
distance between electron and an atom is large com- 
pared to the size of the atom so that exchange interaction 
can be neglected. The diagonalization is important be- 
cause energy differences between quasi-continuum states 
are small. Essentially, all-order methods are necessary, 
at least to include chained self-energy corrections. Pure 
2nd- or 3rd-order expansions will be inaccurate due to 
this reason, and we will illustrate this numerically for 
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TABLE I: "Bag" model artifact. An extra phase shift S due to 
the "bag" boundary condition for an empty cavity of R=15 
a.u. is compared with prediction Qv/ E/2 in the Pauli ap- 
proximation; I is the angular momentum of the state, n is the 
radial quantum number. 



I n. E CCLV 



WE/2 



1 1 4.48 [-2] 1.10 [-3] 1.09 [-3] 

1 3 2.64[-l] 2.69[-3] 2.65[-3] 

1 4 4.39[-l] 3.80[-3] 3.42[-3] 

2 1 7.38 [-2] 1.41 [-3] 1.40 [-3] 
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2nd MBPT, s-wave 



2nd-order MBPT in the next section, but couple-cluster 
methods, which treat some diagrams in all orders, are 
expected to give good accuracy. The simplest future 
improvement for the current BO theory is to take into 
account screening, which is more significant in heavier 
noble-gas atoms. 

Relativistic effects can be also carefully considered, if 
necessary. One effect is the difference in energies be- 
tween for example P1/2 and p 3 / 2 states, which for low- 
energy scattering is small, but becomes more pronounced 
at higher energies. The self-energy correction is also 
slightly different in non-relativistic and relativistic cases 
because intermediate states in the summation are differ- 
ent. One interesting consequence of the use of relativis- 
tic basis in calculations is that the boundary condition 
is not P{R) = 0, but rather P(R) = Q(R), where P(r) 
and Q(r) are large and small components of the radial 
Dirac wavefunction. This is so-call "bag" boundary con- 
dition which is required to avoid Klein paradox [l2j and 
spurious solutions observed in Ref. 0. Using the Pauli 
expansion, it can be shown that the difference in bound- 
ary conditions produces additional phase shift equal to 
a \J E/2, where a is the fine-structure constant and E is 
the energy of the electron. This shift can be obtained if 
we compare energies generated in the empty cavity with 
energies expected from the zeros of the spherical Bcsscl 
functions as illustrated in Table [i] Apparently, the Pauli 
approximation explains well and predicts accurately the 
"bag" shift as long as a^/E/2 < 1. When the shift is 
large, it is necessary to subtract it or even to reanalyze 
this method more carefully, fn the calculations presented 
below the energies were small enough to neglect this ef- 
fect as well as some other relativistic effects. 

The results of our calculations for elastic cross section 
on argon is shown in FigQ] Close agreement with ex- 
perimental data is achieved in the range below 10 eV if 
the self-energy chain corrections are included (BO ener- 
gies are used) and the cross sections from partial waves 
with I — — 2 are added. To emphasize the importance 
of I > contributions in Fig we also plot s-wave cross 
section separately and in Fig|21we compare phase shifts 
from s-, p-, and d-waves. To check that our predictions 
for phase shifts are correct, we compare them with ex- 
perimental phase shifts. Contributions from higher order 
partial waves are much smaller, but can be in principle 
included. At low energies, the dominant contribution 
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FIG. 1: Low energy argon cross section. Comparison of var- 
ious theoretical approximations and of the final accurate BO 
values with experiment. The solid line (partial waves 1=0-2) 
and the dashed line (1=0) show our BO cross sections obtained 
after diagonalization of the Hamiltonian matrix that contains 
self-energy defined by Eq.@. The dotted line is cross section 
obtained from 2nd order MBPT energies without diagonaliza- 
tion. The dash-dotted line shows cross section obtained from 
DHF energies. Experimental results are taken from Refs. [l4| 
(Expt. 1) and (Expt. 2). 



3.0 r 
2.5 - 
2.0 - 




5 10 15 20 

Energy (eV) 



FIG. 2: Comparison of our BO calculations (the solid line) 
with experiment [l6| (the points with error bars) for phase 
shifts of electron scattering from argon. 



comes from s-waves, which is expected; however, at en- 
ergy about 0.36 eV, s-wave sin (5 crosses zero, resulting in 
a minimum of the cross section. In this region the p-wave 
and d-wave contributions become particularly important 
and affect the shape of the Ramsauer-Townsend mini- 
mum. 

Although the diagonalization does not change much 
energies of the quasi-continuum states, the phase shifts 
and cross sections obtained from energies before (dotted 
line) and after (solid line) diagonalization are quite dif- 
ferent, see Fig^ The agreement is achieved only in the 
last case. In the case when argon cross section is cal- 
culated from DHF energies and thus polarization effects 
are ignored, the cross section is completely inaccurate. 
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FIG. 3: Low-energy krypton cross section. The solid line 
shows our theory and dashed line shows pseudopotential cal- 
culations by Plenkiewicz et al. ITsll . Experiment: open circles 
with larger error bars show the cross sections of Gus'kov et al. 
[T^l . solid circles with smaller error bars show cross sections 
of Buckman and Lohmann J20J , and solid circles without error 
bars show results of Subramanian and Kumar 1211 . 



At low energy DHF scattering length is exactly opposite 
to correct value. DHF potential cross section does not 
depend much on energy and this potential in the range 
below 3 eV can be approximated by an infinite poten- 
tial at R < Ro, where i?o =1.42-1.53 a.u. approximately 
equal to the size of the argon electron cloud about 1.56 
a.u. Simple interpretation of this is the repulsion due to 
Pauli exclusion principle. 

Because the experiments at very low energy are dif- 
ficult, we also find scattering length by extrapolating 



our results to zero energy, R sca tt = — 1.47 ± 0.03 a.u. 
There are several other calculations of the argon scatter- 
ing length: -1.63 a.u. by Asaf et al. 17] from studies of 
perturbed optical absorption in gases , -1.492 by Buck- 
man and B. Lohmann ^5| and -1.449 by Ferch et al. [Isj 
from TSC studies, and -1.486 by Saha from low en- 
ergy calculations. Our value disagrees only with the value 
from Ref.[l7j. 

Our theoretical cross section for krypton is shown in 
FigED The agreement with the cross section obtained 
by Plenkiewicz et al. [l9l | from a pseudopotetial is very 
close. The agreement with experiment is also relatively 
good in all range of energies shown, although some dis- 
agreement can be seen near the cross-section minimum, 
which can be due to the inaccuracy of both theory and 
experiment. 

In this paper, we proposed a simple method for cal- 
culations of phase shifts from energies of quasicontin- 
uum states and illustrated its high precision with MBPT 
calculations. The method in general can be applied to 
many scattering problems: electron scattering on various 
atoms and ions, positron scattering, atom-atom scatter- 
ing; however, in each case some specific atomic structure 
method has to be developed to achieve practical preci- 
sion. Discussed BO approximation can be used only for 
electron scattering on closed-shell atoms and ions. 

The author is grateful to Dr. Kuzma for finding rel- 
evant references and to Prof. Happer for discussion of 
physics of electron scattering, for reading the manuscript, 
and for suggestions for its improvement. 
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